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Abstract 

The aim of this paper is to extend the method of jo to scattering problems. We construct 
very effective near-cloaking structures for the scattering problem at a fixed frequency. These 
new structures are, before using the transformation optics, layered structures and are designed 
so that their first scattering coefficients vanish. Inside the cloaking region, any target has 
near-zero scattering cross section for a band of frequencies. We analytically show that our 
new construction significantly enhances the cloaking effect for the Ifelmholtz equation. 
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1 Introduction 

The cloaking problem in electromagnetic wave scattering is to make a target invisible from far- 
field measurements. The difiiculty is to construct permeability and permittivity distributions 
of a cloaking structure such that any target placed inside the structure has zero scattering cross 
section. Since the pioneer works , [14] , and |19j , extensive research has been done on cloaking in 
electromagnetic scattering. We refer to [TU], |^, and pi] for recent development on electromagnetic 
cloaking. One of the main tools to obtain cloaking is to use a change of variables scheme (also 
called transformation optics). The change of variables based cloaking method uses a singular 
transformation to boost the material property so that it makes a cloaking region look like a point to 
outside measurements. However, this transformation induces the singularity of material constants 
in the transversal direction (also in the tangential direction in two dimensions), which causes 
difficulty both in the theory and applications. To overcome this weakness, so called 'near cloaking' 
is naturally considered, which is a regularization or an approximation of singular cloaking. In jl3| . 
instead of the singular transformation, the authors use a regular one to push forward the material 
constant in the conductivity equation, in which a small ball is blown up to the cloaking region. 
In [12] , this regularization point of view is adopted for the Helmholtz equation. See also [HI [18] . 
It is worth mentioning that there is yet another kind of cloaking in which the cloaking region is 
outside the cloaking device, for instance, anomalous localized resonance [16l[l7l|3]. See also [2]. 

The purpose of this paper is to propose a new cancellation technique in order to achieve 
enhanced invisibility from scattering cross section. 
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Our approach extends to scattering problems the method first provided in [(3 to achieve near- 
cloaking for the conductivity problem from Dirichlet-to-Neumann measurements. It is based on 
the multi-coating which cancels the scattering coefficients of the cloak. We first design a structure 
coated around a perfect insulator to have vanishing scattering coefficients of lower orders and show 
that the order of the scattering cross section of a small perfect insulator can be reduced significantly. 
We then obtain near-cloaking structure by pushing forward the multi-coated structure around a 
small object via the standard blow-up transformation. Note that such structure achieves near- 
cloaking for a band of frequencies. In order to illustrate the viability of our method, we give 
numerical values for the permittivity and permeability parameters and the thickness of the layers 
of scattering coefficients vanishing structures. 

This paper is organized as follows. In the next section we derive the multi-polar expansion 
of the solution to the Helmholtz equation, and define the scattering coefficients. In Section 3 we 
characterize the scattering coefficients vanishing structures. In Section 4 we show that the near- 
cloaking is enhanced if a scattering coefficients vanishing structure is used. In Section 5 we present 
some numerical examples of the scattering coefficients vanishing structures. 

Even though we consider only two dimensional Helmholtz equation in this paper, the same 
argument can be applied to the equation in three dimensions. The multi-coating technique devel- 
oped in this paper can be extended to the full Maxwell equations, which will be a subject of a 
forthcoming paper. 

2 Scattering coefficients 

For fc > 0, the fundamental solution to the Helmholtz operator (A-f fc^) in two dimensions satisfying 

(A-f fc2)r,(x) = 5o 

with the outgoing radiation condition is given by 

r,(x)^-^i/«(fc|x|), 

where H^^^ is the Hankel function of the first kind of order zero. For |x| > |y|, we recall Graf's 
addition formula [22], 

i7^^)(fc|x -y\)^Y. Hi^Hkl^De-'-'-JniklyDe-^'y, (2.1) 

where x = (|x|, 0x) and y — (|y|, 9y) in polar coordinates. Here Hn ^ is the Hankel function of the 
first kind of order n and J„ is the Bessel function of order n. 

We first define the scattering coefficients of an inclusion (with two phase electromagnetic ma- 
terials) and derive some important properties of them. 

Let D he a, bounded domain in with Lipschitz boundary dD, and let (eo,/io) be the pair of 
electromagnetic parameters (permittivity and permeability) of \ -D and (ei,/xi) be that of D. 
Then the permittivity and permeability distributions are given by 

e = eoxi^^ \ D) + eix{D) and fi ^ fioxi^'' \ D) + fnxiD) . (2-2) 

Here and throughout this paper x{D) is the characteristic function of D. 

Given a frequency w, set k = w^/ei/TT and fco = w^eo/io- For a function U satisfying (A-t-fep)?/ = 
in M^, we consider the scattered wave u, i.e., the solution to the following equation: 

V-ivu + w^eu = inM^ 

M (2.3) 

(u — U) satisfies the outgoing radiation condition. 
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Let 5^[<(5] be the single layer potential defined by the kernel Ffc, i.e., 

5^M(x)= / rfe(x-yV(y)da(y). 

JdD 



IdD 

Let S'^ be the single layer potential associated with the kernel F^^ . Then, from, for example, [5] , 
we know that the solution to (|2.3p can be represented using the single layer potentials 5^° and iS|, 
as follows 



^ fc/(x)+4"[V.](x), i,eR^\D, 

where the pair (t/J, V') G L'^{dD) x L'^[dD) is the unique solution to 



(2.4) 



ld{S%M) 



fi dv 



1 d{s^M) 



1 a[/ on5i:i. (2.5) 

Mo 



Here, + and — in the subscripts respectively indicate the limit from outside D and inside D to dD 
along the normal direction and d/dv denotes the normal derivative. It is proved in [5] that there 
exists a constant C ~ C{k, ko, D) such that 

M\maD) + \\M\L-[dD) < C{\\U\\L2(^gD) + |1 VC/|U2(gc))- (2.6) 

It is also proved in the same paper that there are constants po and C = C{k,kQ, D) independent 
of p as long as p < po such that 

W^ph^idD) + Uph^idD) < C{\\U\\mOD) + ||Vf/||L2(gc)), (2.7) 

where {ipp,ipp) is the solution of (|2.5p with k and /cq respectively replaced by pk and pk^. 

Note that the following asymptotic formula holds as |x| — >■ oo, which can be seen from ()2.4p 
and Graf's formula (12.11): 



u{^) - L/(x) - E i/(i)(fco|x|)e-^=' / J„(fco|y|)e-™^^^(y)da(y). (2.8) 

Let [iprm'^m) be the solution to (12. 5p with Jm(fco|x|)e™^^ in the place of J7(x). We define the 
scattering coefficient as follows. 

Definition 1. The scattering coefficients Wnm, m,n (z Z, associated with the permittivity and 
permeability distributions e,/i and the frequency uj (or k,kQ,D) are defined by 



Wnm = Wnm[e,ti,u;]:= (^0 Iy I )e^"''- V'm (y)da(y) . (2.9) 

JdD 

We first obtain the following lemma for the size of 
Lemma 2.1. There is a constant C depending on {e,p,uj) such that 

\Wnm[e,li,u}]\ < for alln^meZ. (2.10) 

Moreover, there exists po such that, for all p < po, 

fj\n\ + \m\ 

|W^„™[e,M,HI < foralln,meZ, (2.11) 

where the constant C depends on (e, p, uj) but is independent of p. 
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Proof. Let C/(x) = Jm(fco|x|)e™^^ and {(pm,'ipm) be the solution to (I2.5p . Since 

'^'"W^^^^O^ 2.12 
\/27r|m| V2|m|/ 

as m — C!0 (see [1]), we have 



|to|I 

for some constant C. Thus it follows from (12.61) that 



UmWmaD) < (2.13) 

for another constant C. So we get (|2.10p from (|2.9p . 

Now let ((/?,„, V'm) be the solution to (|2.5I) with fc and fcp replaced by pfc and pkQ. One can see 
from ()2.7p that (|2.13p still holds for some C independent of p as long as p < po for some po- Note 
that 



Wn„,hp,puj]^ / J„(pfco|y|)e ^^^V'm,p(y)rfcr(y), 

where {(pm,p, 4'm,p) is the solution to (|2.5p with fc and /co respectively replaced by pk and pfco and 
^m(fcop|x|)e™^'' in the place of U{x). So one can use (|2.12p to obtain ()2.1ip . This completes the 
proof. □ 

If U is given as 

C/(x) = am(C/)Jrn(fco|x|)e™^= (2.14) 

where am(U) are constants, it follows from the principle of superposition that the solution {(p,ip) 
to (|2.5p is given by 

mGZ 



Then one can see from (|2.8p that the solution u to (|2.3p can be represented as 

7/(x) - [/(x) = ^ i/,W(fco|x|)e"''= ^ W^„„a™(C/) as |x| ^ oo. (2.15) 

In particular, if U is given by a plane wave e***^"^ with k • k = /cq, then 

zi(x)-e^'^-'' = -^^i/«(A;o|x|)e'"''='^PF„„e^"(^-''>') as |x| ^ cx), (2.16) 

where k = fco(cos0k,sin6'k) and x = (|x|,(?x). In fact, since 

g»feorsi„e ^ ^ J,„(fcor)e™^ 

we have 

gik-x ^ ^ e""(^-''>'V„(fco|x|)e''"*^ (2.17) 

and 

= E (2.18) 

mGZ 
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Thus (I2.16|) holds. It is worth emphasizing that the expansion formula (|2.15p or (I2.16P determines 
uniquely the scattering coefficients Wnm, for n,m G Z. 

We now show that the scattering coefficients are basically the Fourier coefficients of the far-field 
pattern (the scattering amplitude) which is 27r-periodic function in two dimensions. The far-field 
pattern A^oi^, iJ,,io], when the incident field is given by e''^'^, is defined to be 



it(x)-e'''"^ = -ie T^:^=^Aoo[e,/i,H(6'k,6'x)+o(|x| as |x| c5o. (2.19) 
V87rfco|x| 

Recall that 

^o^^W ^ y'^e*^*-^) ast^oo, (2.20) 

where ^ indicates that the difference between the right-hand and left-hand side is 0{t^^). If |x| 
is large while |y| is bounded, then we have 

|x-y| = |x|~|y|cos(0.-0y) + O(^), 

and hence 



7I^^Vfco|x-y|) - e-^r /^^e'''^"(l^l-|yl™^(^--^^» as Ixl ^ oo. 

V TTfco x| 



Thus, from we get 

w(x) - e**"^ - ^ie-'r-^^^^^^ f e-'ko\y\cos{e^~ey)^(^y^ ^^^^^ |^| ^ 21) 
■\/87rfco|x| JdD 

and infer that the far-field pattern is given by 

A^{e^,e^)^ [ e-'''<'\^\^°<'>-''y^{y) do{y), (2.22) 

JdD 

where -0 is given by (I2.18p . 
Let 

Aoo(^k,^x) = ^fo„(^k)e™^'' 
be the Fourier series of Aoo(^k, •)• From (|2.22p it follows that 

bn{e^.) = ^fj I ^-^'=''1^1 ™-(«-»y)0(y) da{y) e-"''- d0^ 



I 

2^ 



dD Jo 



dD 
2-K 



e 



Since ^ 

_L f e-*'=''lyl'=°"(''--''^)e-"''-d6lx = J„(fco|y|)e~™(''^+*\ 
27r 







we deduce that 



bn^^ I J„(fco|y|)e-™<^-+5V(y)rf^(^y)- 

IdD 



Using (|2.18p we now arrive at the following theorem. 

Theorem 2.2. Let 6 and 9' be respectively the incident and scattered direction. Then we have 

A^[e,^l,uJ]{e,e')= ^ i("~")e"'«V„„[e,/.,c.]e-*™^ (2.23) 
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We emphasize that the series in (|2.23|) converges uniformly thanks to (|2.10p . We also note that 
if U is given by (|2.14p then the scattering amplitude, which we denote by Aoo[£, fJ-,ij-']{U,9'), is 
given by 

A^[e,fi,uj]{U,e') = ^i-'^e™"' ^ iy„™a,„(C/). (2.24) 

The conversion of the far- field to the near field is achieved via formula (|2.16l) . 
The scattering cross section 5[e,/i,aj] is defined by 



dO. 



(2.25) 



See [3 [20] . As an immediate consequence of Theorem 12.21 we obtain the following corollary. 
Corollary 2.3. We have 



S[e, fx, w] (0') = 271 ^ J2 i'^'Wnm [e, w]e"«' 



(2.26) 



It is worth mentioning that the optical theorem [71 [50] leads to a natural constraint on Wm. 
In fact, we have 



^mAoo[e,fi,u;]{0',9') = -J — S[e,fi,u]{e'), V 0' e [0,2^], 



(2.27) 



or equivalently, 



(2.28) 

In the next section, we compute the scattering coefficients of multiply coated inclusions and 
provide structures whose scattering coefficients vanish. Such structures will be used in Section 
[4] to enhance near cloaking. Any target placed inside such structures will have nearly vanishing 
scattering cross section S, uniformly in the direction 9' . 



3 S-vanishing structures 

The purpose of this section is to construct multiply layered structures whose scattering coefficients 
vanish. We call such structures S-vanishing structures. We design a multi-coating around an insu- 
lated inclusion D, for which the scattering coefficients vanish. The computations of the scattering 
coefficients of multi-layered structures (with multiple phase electromagnetic materials) follow in 
exactly the same way as in Section [5] The system of two equations (|2.5p should be replaced by a 
system of 2x the number of phase interfaces ( — 1 if the core is perfectly insulating). 
For positive numbers ri, . . . , r^+i with 2 — ri > r2 > ■ ■ ■ > J'l+i = 1, let 

yl, := {x : r,+i < < r,}, j = l,...,L, Ao:=R^\A^., Al+i{= D) '.^ {x : \x\ < 1}. 

Let (pj, tj) be the pair of permeability and permittivity of Aj for j=0,l,...,L + l. Set /io = 1 
and eo = 1. Let 

L+l L+l 

^ = X! ^jX(^i) ^'^'^ ^ = ^jX(^i)- (3-1) 

3=0 j=0 
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In this case the scattering coefficient Wnm = W„m[/z, e,aj] can be defined using (|2.15|) . In fact, 
if u is the solution to 

1, 



V • -Wu + oj^eu = 



(3.2) 



with the outgoing radiation condition on u — U where U is given by (j2.14p , then u—U admits the 
asymptotic expansion (I2.15|) with fco = w^eo/io- 

Exactly like the conductivity case jBj one can show using symmetry that 



if TO 7^ n. 



(3.3) 



Let us define Wn by 

Wn Wn„. (3.4) 

Our purpose is to design, given N and w, fi and e so that e,uj] — for |n| < N. We call such 

a structure (//, e) an S-vanishing structure of order N at frequency uj. Since 

H^-l - (-l)"i^«^' and 

J-„ = (— 1)"J„, we have 



and hence it suffices to consider Wnn only for n>Q. 
Note that ((2?28)) leads to 



9to ^ W„[e,Ai,a;] = - J^Xl 



W„[e,^, 



(3.5) 



(3.6) 



Let fcj :— Lj^^jEj for j = 0, 1, . . . , L. We assume that — oo, which amounts to that the 

solution satisfies the zero Neumann condition on |x| = ri+i(= 1). To compute Wn for n > 0, we 
look for solutions it„ to (|3.2p of the form 

(3.7) 
(3.8) 



u„(x) = af^ Jn{kjr)e"'' + Hi^\k,r)e"^\ x G A,, j - 0, . . . ,i, 



with On"'' = 1. Note that 



The solution u„ satisfies the transmission conditions 



= Un\- and 



1 dUn 



dv 



1 dUn 



+ H 



on X = To 



for J = 1, . . . , L, which reads 



Jnjkjrj) 

Jn(kj—irj ) 



(ilW)'(fc,r,) 



J„(fcj_irj) 



















i- 









Mi-i V /^i-i 

The Neumann condition ^^|+ = on |x| — r^+i amounts to 



j;(fci) 

Combining p.9p and (13.101) . we obtain 



(3.9) 











0" 














P(")[6,A*,H 



J")' 

^0 

r,(») 



(3.10) 



(3.11) 
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where 



(«) («) 

1^21 P22 . 



- n 



L 



H 



JnikjTj) 



''-Jj^{kjr^ 





Jn{kj-irj) 
-J'nih-irj) 



In order to have a structure whose scattering coefficients Wn vanishes up to the A^-order, we 



need to have feg"'' = (when Oq"'' = 1) for n = 0, . . . , A^, which amounts to 



p^"^ = for n = 0, . . . , 



(3.12) 



because of p. lip . We emphasize that P22 7^ 0. In fact, if = 0, then (I3.1ip can be fulfihed 



(") _ 



with a^"^ = and = 1. It means that there exists (/u, e) on \ D such that the following 
problem has a solution: 

1, 



V ■ -Vu + uj^eu = mR^\D, 

du „^ 
= on dD, 

+ 



dv 



(3.13) 



U(x) = M'^(fcor)e 



for Ixl = r > 2, 



which is not possible. 

We note that (|3.12p is a set of conditions on {fij,ej) and rj for j = 1, . . . , L. In fact, pji^ is 
a nonlinear algebraic function of fij, Sj and Vj, j = 1, . . . , L. We are not able to show existence 
of {fj,j,ej) and Vj, j — 1, . . . ,L, satisfying (I3.12p even if it is quite important to do so. But the 
solutions (at fixed frequency) can be computed numerically in the same way as in the conductivity 
case [6]. 

We now consider the S- vanishing structure for all (low) frequencies. Let uj be fixed and we look 
for a structure (/i, e) such that 

Wn[n, e, puj] = for aU \n\ < N smd p < po (3.14) 

for some po- Such a structure may not exist. So instead we look for a structure such that 

VK„[/i,e,H = o(p2^) for all |n| < iV and p ^ 0. (3.15) 

We call such a structure an S-vanishing structure of order N at low frequencies. 

To investigate the behavior of Wn [p, e, pto] as p ^ 0, let us recall the behavior of Bessel functions 
for small arguments. As i — > 0, it is known that 



1 



a'' 



Y^{t) 



2"Vr(n + l) r(ra + 2) 2!r(n + 3) 3!r(n + 4) 
(it)-" '^^^ {n-l-\)\A 



E 

i=0 



{-/Y + -H-t)Ut) 



i^y(V(/ + l)+^(n + / + l))> 4^ ' 



1=0 
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where "0(1) = —7 and ■(/'(n) = —7 + X^ILi^ ^'^'^ n > 2 with 7 being the Euler constant. In 
particular, if n = 0, we have 



Ut)^i--e + l-/ + o{t^), 



(3.16) 



Yo{t) = -Ini + -(7-ln2) - —t^lnt+i- ^-{-f -\TL2)\t^ + Oif^lnt). (3.17) 



Plugging these formulas into (|3.1ip . we have 



xn 











-YP + 0(p3) J^p-^ + 0(plnp) 



p i+0(plnp) — 

TTCUfljrj TT^ 



1 



1 \ In p 



(jjfij-irj ^IJ-jTj J p 



+ o{p-') 




2i TT- 



P + 0(p3) 



0(p) 



2i 



TTOJpj-irj 



-p ^ + 0{p\np) 



(3.18) 



and, for n > 1, 



xn 



i=i 

^z2"r(,i + l) „ 1 











' , p"-^+0(p") »2"r(n + l) 1 
2"r(n + l) ^ 7rfc"+i ^ VP ; 



„-7^+l^ 



Pj 2"r(n + 1)' 



7r(fcjrj)" 



0(p") 



2"r(n + 1) 



1 



2"r(n + l) 




p" + 0(p"+') 

p 



2"r(n + 1) 
i2"r(n) 



+0(p") 



p" + 0(p"+i) 

p - ^ + o(p-") 



7r(fcj_irj) 
T-^ z2»r(n + 1) ^_„_, 
Mi-i 7r(fcj_irj)"+i 





xR 



_2"r(n + 1)"^ 

a,(6, + 1) + 0(1) c,{h, - l)p-2" + o(p-2"y 

fe,+l 



where 



0(1) 



(3.19) 



, 6j 



Pj _ i22"r(n)r(n + 1) 



with kj = ujyfejpj. 

From the above calculations of the leading order terms of [e, p, pui] and the expansion 
formula of J„(t) and Yn{t), we see that P21'* cind P22'* admit the following expansions: 



{N-n) L + 1 

p^-\p,e,t)=e~'\f^-\p,e)+ E4?(a^'^)^''(1"^)''+«(^'''"'") 

1=1 j=Q 



(3.20) 
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and 



{N-n) L+1 



2N-2n\ 



1=1 j=0 



for t ^ pcu and some functions /o"\5o"^ //i^' ^^'^ d'fj' independent of t. 
Lemma 3.1. For any pair of e), we have 



(3.21) 



(3.22) 



PtooL For n = 0, it follows from ((3?T8| that 



(0)/ N 2i 



n 



^0. 



Suppose n > 0. Assume that there exists a pair of {p,e) such that g^\^,e) = 0. Then the 
solution given by p.7p with ap"-* = and ^q""* = 1 satisfies 



V • -Vu + p^w^eu = in \ D, 
M 
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(3.23) 



[u(x) = Hi^\pkor)e''^'' for |x| = r > 2. 



Let w(x) := limp_^o p"'u(x). Then v satisfies 



V • -\/v = in R2 \ D, 
= ondD, 



dv 
i;(x 



(3.24) 



i2'T(n) 



■r-"e™^ for Ixl = r > 2, 



□ 



^(eoAio)"/2 

which is impossible. Thus gj^\fi, e) 7^ 0, as desired and the proof is complete. 

Equations p.20p and p.2ip together with the above lemma give us the following proposition. 
Proposition 3.2. For n> 1, let Wn be defined by p^P- We have 



I (N-n) M^j \ 

Wn[p,e,t]{0,0')^t'" W^[p,e]{e,e')+ Y W^^[p,e]{e,e')t^'{lnty +o(i2^) (3.25) 

\ 1=1 j=0 / 

uniformly in (6, 9'), where t = puj, Mn^i : {L + l)(iV — n) (L being the number of layers), and the 
coefficients e] and Wli^[p,€] are independent oft. 

To construct an S- vanishing structure of order N at low frequencies, we need to have a pair 
(/i, e) of the form p.ip satisfying 

W°[p, e] = 0, and Wl^^[p, e] = for < n < iV, 1 < ^ < (A^ - n), 1 < j < M„,,. (3.26) 

We construct numerically such structures for small N in the last section. 
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4 Enhancement of near cloaking 

In this section we show that the S- vanishing structures (after a transformation optics) enhance the 
near cloaking. 

Let e) be a S- vanishing structure of order N at low frequencies, i.e., p.26p holds, and it is 
of the form It follows from ((2lT1) . Theorem [l^ and Proposition that 



uniformly in (9, 6') ii p < po for some po. 
Let 



1 



*i(x)==-x, xeM^ 

p 



Then we have 



p o 1, e o vl/j. , u! 



Aoc[p, e, puj]. 



(4.1) 

(4.2) 
(4.3) 



To see this, let u be the solution to 

1 



V • 



(AiO*i)(x) 



Vu(x) +a;2(eo ^'i)(x)u(x) = in \ Sp, 



— = on dBp, 
ou 



(4.4) 



, {u ~ U) satisfies the outgoing radiation condition, 
where t/(x) = g»feo(cose,sine) x^ jj^^^.^^ 5^ is the disk of radius p centered at 0. Define for y = ^x 



i2(y) :=( uo*7M(y) =Uo^'J(y) and (/(y) - ( (7 o J (y) 



Then, we have 



V • 



M(y) 



Vyt2(y) + p'w'e(y)?i(y) = h 



^ = onaSi, 

{u — U) satisfies the outgoing radiation condition. 
From the definition of the far-field pattern A^o, we get 



(4.5) 



^. gifco|x| 

{u - [/)(x) ie~T_==__A^ 



V87rfco|x| 



p o vj/j^ ^ £ o ^! 1 , w 



?') as |x| — T' oo. 



and 



gipfcolyl 

»e~T^— Aoo[m, e,pi:j](6', g^) as |y| ^ oo. 



\/87rpfco|y 

where x = |x|(cos0',sin6''). So, we have (|4.3p . It then follows from (|4.ip that 



/i o vpj, , e o ^ J. , w 



9') = o{p^^) 



We also obtain from (|2.26p 



S 



(4.6) 
(4.7) 



It is worth emphasizing that (/iovl/j_^£ovl/j_)isa multi-coated structure of radius 2p. 

We now apply a transformation to the structure {po"^ i,to Vfj.) to blow up the small disk of 
radius p. Let us recall the following well-known lemma (see, for instance, [TO])- 
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Lemma 4.1. Let F be a diffeomorphism o/ onto such that -F(x) is identity for |x| large 
enough. If v is a solution to 

V ■ A\7v + uj^qv = mR^, 
subject to the outgoing radiation condition, then w defined by w(y) = v{F^^(y)) satisfies 

\7 ■{F^A)S/w + uj^{F^q)w:^0 m M^ 

together with the outgoing radiation condition, where 

DF{x)A{x)DF^{x) , , , , , q(x) 



(4.8) 



and {F„q){y) 



det DF{x) ' det DF{x) 

with X = F^^{y) and T being the transpose. 

For a small number p, let Fp be the diffeomorphism defined by 

for |x| > 2, 



(4.9) 



^p(x) 



3-4p 



1 



2(1 -p) 4(1 -p)' 
1 
2 



for 2p < |x| < 2, 
for p < |x| < 2p, 
for |x| < p. 



(4.10) 



We then get from (|4.6p . (I4.7p . and Lemma [4. II the main result of this paper. 

Theorem 4.2. // (/i, e) is a S-vanishing structure of order N at low frequencies, then there exists 
Pa such that 



and 



S 



{Fp),ipo^^),{Fp),ieo^^),Lo {9,9')^o{p'^) 



{FpU^iO^>^),{FpU^O^^),UJ (0') = o(p4^), 



(4.11) 
(4.12) 



for all p < po, uniformly in 9 and 9' . 



We first note the cloaking enhancement is achieved for all the frequencies smaller than uj. This 
is because (|3.15l) holds if we replace w by w' < cj. Then it is worth comparing (I4.1ip with (|4.ip . In 
()4.ip . (/i, e) is a multiply layered structure between radius 1 and 2 in which each layer is filled with 
an isotropic material, and enhanced near cloaking is achieved for low frequencies puj with p < po- 
On the other hand, in (I4.11[) the frequency lu does not have to be small. In fact, (14. lip says that 
for any frequency uj there is a radius p which yields the enhanced near cloaking up to o{p'^^). 
We emphasize that (-Fp)* (/i o ) and (F„)^{e o vj/ 1 ) are anisotropic permittivity and permeability 
distributions. It is not clear whether we can achieve enhanced near cloaking at high frequencies 
by using isotropic layers as done in (|4.ip . 



5 Numerical examples 

In this section we provide numerical examples of S-vanishing structures of order N at low frequen- 
cies, i.e., structures (//, e) of the form p.ip satisfying p.26p . To do so, we use the gradient descent 
method to minimize over (//, e) the quantity 



(N-n) M„,, 



1=1 j=0 



(5.1) 
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where W^°[^, e] and W^^^fi,e] are the coefficients of W„[/i, e,i] in p.25p . These coefficients are 
expressed in terms of Bessel functions and their derivatives. We use the symbohc toolbox of 
MATLAB in order to extract these coefficients. 

It is quite challenging numerically to minimize the quantity in (jS.ip for large N. In numerical 
examples we take iV = 2. In this case one can show through tedious computations that the nonzero 
leading coefficients of M^„[/i, e, t] are as follows: 

• [t^,t'^,t'^logt] for n = 0, I.e., W^'", Wo'°, W^'^; 

• for n = 2, i.e., . 

Results of computations are given in Figure [01 when we use 0, 1, 2 layers {L — 0,1,2). 
The computed material parameters are fix — 0.6 and e = 4/3 when we use one layer, and /i = 
(1.4905,0.27594), e = (1.09271, 1.6702) when we use two layers. The first two columns show the 
distribution of fi and e for L = 0, 1, 2. The last column shows the values of coefficients of W„[^, e,t]: 
Wo'°, Wo'°, indicated by (0,0,0); W^, W^'", Wl'^ indicated by (1, 1, 1); indicated by 2. 

The figures show that if L = then none of the coefficients is zero; \i L = 1, then the coefficient 
of is close to zero; ii L ^ 2 all the coefficients are close to zero. 

We then compute the scattering coefficients e, t] using the computed {fi, e) for L = 0, 1, 2. 

Figure [521 shows the results of computations for n = 0, . . . , 4 and t = 1, 0.1, 0.01. It clearly shows 
that Wn[n,€,t] for n < 2 gets smaller as the number of layers increases. 

6 Conclusion 

We have shown near-cloaking examples for the Helmholtz equation. We have designed a cloaking 
device that achieves enhanced cloaking effect. Any target placed inside the cloaking device has an 
approximately zero scattering cross section. Such cloaking device is obtained by the blow up using 
the transformation optics of a multi-coated insulating inclusion. In the numerical example, to 
minimize the scattering coefficients of the multi-coated insulating inclusion up to the second order, 
we set 2 layers with different permittivity and permeability properties. The technique presented in 
this paper can be extended to full Maxwell's equations. This will be the subject of a forthcoming 
paper. 
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